Exact solution of a Brownian inchworm model for self-propulsion 



o 
o 



I 



o 



X 



A. BauleE 

The Rockefeller University, 1230 York Avenue, New York, NY 10065, USA 

K. Vijay Kumai0 and Sriram RamaswamjH 

CCMT, Department of Physics, Indian Institute of Science, Bangalore 560012, India 

We present the exact solution of a Brownian inchworm model of a self-propelled elastic dimer 
which has recently been proposed in [K. V. Kumar et al, Phys. Rev. E 77, 020102(R) (2008)] as 
a unifying model for the propulsion mechanisms of DNA helicase, polar rods on a vibrated surface, 
crawling keratocytes, and Myosin VI. 
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I. INTRODUCTION 



^ , Methods of nonequilibrium statistical mechanics are particularly useful in describing transport processes in biological 
^ ■ systems. In this article we investigate a nonequilibrium model of a molecular motor (or self-propelled particle), in 
which unidirectional movement is generated by rectification of Brownian diffusive motion. A variety of dynamical 
processes in cells arise from motor proteins such as kinesin, myosin or DNA helicase, which are driven by hydrolysis 
of ATP into ADP and move along tracks consisting of microtubules, actin filaments, or DNA respectively [U, 0, [Sj]- In 
. larger assemblies molecular motors give rise to motion of bacteria and are ultimately responsible for muscle contraction 
and macroscopic movement of living organisms. But such directed motion is also observable in the non-living world, 
for example in agitated granular matter 0, Q. From a physicist's point of view, one of the main interests in these 
C ' self-propelled systems concerns the basic principles by which directed motion is generated from isotropic energy input 

[ in the absence of an externally imposed gradient. 
(-H I In the above mentioned systems the unidirectional movement of the center of mass results quite generically from a 
Q ' coupling to the internal coordinates under the influence of (i) energy input due to chemical or mechanical nonequi- 
O . librium noise and (ii) an asymmetric environment for the internal coordinates. A widely used model of a Brownian 
motor implementing these principles is based on an asymmetric ratchet potential, acting on the centre of mass of 
, the particle Models of this type date back to pioneering work of Smoluchowski, Feynman, and Huxley (see 0,[1| 
J> ' and references therein), and have subsequently attracted a lot of research activities. The basic concept is that of 
! a Brownian particle moving in an asymmetric periodic potential. Due to the second law of thermodynamics, the 
asymmetric energy landscape and thermal noise alone are not sufficient to generate directed motion. A common way 
to model the effect of the required nonequilibrium energy input is to assume that it gives rise to a rapid increase in 
. - the system temperature thus effectively resulting in a periodic switching between high and low temperature states. 
CO At high temperature the particle can freely diffuse over the potential barriers, whereas at low temperature it will be 
, trapped in a potential well. Due to the spatial asymmetry a periodic switching between these states will cause a net 
movement in a direction prescribed by the asymmetry. Such ratchet models might serve as crude simplifications of 
the movement mechanism of motor proteins, where the asymmetric periodic potential results from the interaction of 
the motor with the track and ATP hydrolysis is the cause of the sudden temperature jumps. 

An alternative picture of a Brownian motor can be established by modelling the motor-track interaction in terms 
of an effective friction force 3. Very recently an inchworm model of a self-propelled particle has been proposed 
' in [l^, where the required spatial asymmetry is implemented by such a friction force. This model consists of two 
beads connected by an elastic spring, which are driven by thermal and nonequilibrium noise. Each bead experiences 
a different, stretch dependent friction, which is sufficient to generate a non-zero center of mass velocity in a fixed 
direction. On the basis of this new rectification mechanism the propulsion in systems as diverse as polar rods on a 
vibrating surface, DNA helicase, and Myosin VI on actin, can be understood in a unified way, as discussed in p^. 

The purpose of the present paper is to further investigate the new Brownian inchworm model of Ref. [lo| . which 
so far has only been treated perturbatively and numerically. In particular, we present exact solutions for the average 
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FIG. 1: Depiction of the Brownian inchworm. A movement as shown here can be observed, for example, if (i) the friction on 
bead 2 is constant and stretch-dependent on bead 1 (increases for larger stretch), and (ii) the nonequilibrium noise-strengths 
are equal. 



center of mass velocity and the distribution of internal coordinate of the inchworm in the overdamped regime. The 
exact solution allows for a strict derivation of some of the perturbative results and shows excellent agreement with 
simulation data. Furthermore, the effect of a more realistic spring potential, such as the FENE spring, is investigated. 



II. THE BROWNIAN INCHWORM 



The Brownian inchworm is described in terms of two coupled Langevin equations for the positions of the two beads 
xi and X2 under the influence of thermal and nonequilibrium noise [lOi] : 

m,i,{t) + -f^ix)x^{t) = -~d,U{x) + ^j2-1,{x)kBTlf^i{t) + ^/MQ^{t). (1) 

Here, di = d /dxi with i G {1,2} and U (x) is the spring potential depending on the relative coordinate x = xi—X2 — Xf) 
for an equilibrium extension xq. The crucial ingredient in this model is the stretch dependent friction 7i(x) > 0, acting 
independently on the two beads. The thermal noise 7]i{t) has Gaussian characteristics with zero mean {rii{t)) = 
and correlation {'rii{t)rjj{t')) — 5ij5{t — t'). The thermal character is expressed in the prefactor y^2'-fi{x)kBT of 
rii{t), stemming from the fluctuation dissipation relation. Likewise, we assume Gaussian statistics with zero mean 
and delta-correlation for the nonequilibrium noise which, in contrast to the thermal noise, does not satisfy a 

fluctuation-dissipation relation as it is considered to represent the external nonequilibrium energy input. 

The two uncorrelated Gaussian processes in Eq. ([!]) can be superposed into a single one, denoted by ^i{t), whereby 
— s-iid {£,i{t)£,j{t')) — SijS{t — t'). The equations of motion then read 

miXi{t) +^i{x)xi{t) = ~diU{x) + gi{x)£,i{t), (2) 

where 

g^ix) = ^2j,{x)kBT + A,. (3) 

Important time scales in this system are the relaxation time of the center of mass velocity and of the fast and 
slow modes of the relative coordinate. When the center of mass relaxes much faster than the relative coordinate, 
the system is effectively in the overdamped regime and the influence of inertia can be neglected. Independent of the 
precise values of the timescales, it is clear that, in a system where dissipation dominates, the correlation time of the 
noise is the shortest, that of the velocities is next ni/j), and that of the relative coordinate is the longest (~ 7/k) 
where k is the spring constant, 7 is a typical value of the x-dependent damping, and m is a typical mass. Thus one 
is always in the regime where the adiabatic elimination procedure can be done with a white noise approximation. 

In the steady-state regime the Brownian inchworm exhibits unidirectional movement which can be described as 
follows. Let us consider the case of constant friction on bead 2 while bead 1 experiences stretch-dependent damping 
(increasing for larger stretch). Gaussian noise results in a high frequency of stochastic kicks compared with the 
relaxation rate of the relative coordinate, such that the dimer basically never relaxes while being under the influence 
of the noise. Since bead 2 has a higher mobility it will accumulate a larger displacement than bead 1 whose movement 
is inhibited with increased stretch. Ultimately bead 1 follows bead 2 since this direction is favoured due to the friction 



3 



and the spring force. In this case the net effect of the fluctuations is a movement of the center of mass in the negative 
direction (see Fig.[T|). 

However, both center of mass velocity and direction of movement depend on the magnitudes of the noise strengths 
Ai and A2. For A2 > Ai a current reversal can be observed. We find quite generally that, whenever A2 ^ ^1, the 
direction of movement is dominated by the asymmetric noise input rather than by the asymmetric friction. This is 
further discussed in Sec. IV A] below. 



III. EXACT SOLUTION IN THE OVERDAMPED REGIME 



An analytical treatment is accessible in the overdamped regime where the dynamics exhibits a time scale separation 
between position (slow) and velocity (fast) degrees of freedom. It is important to note that the position-dependent 
friction in the equations of motion ^ requires a careful inspection of the adiabatic elimination See [13, [3| 

for a discussion of the issues involved. At least for white noise sources, the mathematically sound way to carry 
out the elimination of the velocity is to start with the Fokker-Planck equation for the distribution of position and 
velocity and extract the Smoluchowski equation for the distribution for position alone, through an expansion in I/7. 
Working backwards from the resulting Smoluchowski equation yields the overdamped Langevin equation [ij, El] for 
the position coordinates Xi (see Appendix E| 



7i(a;) 2 li{xY -fi{x) 

with multiplicative noise interpreted a la Stratonovich. 

We will see below that Eq. ([4]) with Stratonovich interpretation gives rise to the correct equilibrium distribution 
for the relative coordinate x, namely Pst{x) the absence of the active noise. If one wants to attribute 

a different interpretation to the multiplicative noise terms, an additional drift term has to be added in Eq. Q in 
order to obtain the physically correct distribution for pst{x). There exists in fact a whole one parameter family of 
overdamped Langevin equations which are equivalent to Eq. ^ . 

In the following the prime denotes a derivative with respect to the relative coordinate x, so that diU{x) — 
—d2U{x) = U{x)' . With Eq. ([4]) the equation of motion for the center of mass coordinate can be written as 



2-fi{xy 272 (x) 

Likewise, the equation of motion for the relative coordinate reads 



2 V7i(a^) l2{x) J 4 V 71 (a;)^ 72(2;)^ 



71 (x) l2{x) J 2 V li{xY 72(2;) 

+^a(^)-^6W. (6) 

71 (x) 72 (x) 

The key mathematical observations here are: (i) the equation of motion for the center of mass is a function of the 
relative coordinate and the noise only, (ii) The equation of motion for the relative coordinate is independent of the 
center of mass coordinate. As a consequence of these properties we obtain the average center of mass velocity directly 
by averaging Eq. ([5|) over the stochastic realizations 



2 \V7i(a;) 72(2;)/ / 4\V 7i(a;)^ l2{xY 

'^''^%it)) + U'-^Ut)). (7) 



2 \7i(a;) 7 2 \72(x) 

It is important to note that the contributions of the multiplicative noise terms are non-zero due to the Stratonovich 
interpretation and can be calculated with the help of Novikov's theorem for Gaussian processes. 

For an arbitrary functional u[?7(t)] of a delta-correlated zero mean Gaussian process r]{t), this theorem states that 
the average over the product {u[ri{t)]ri{t)) can be expressed as an average over the response of u to a change in r]{t) 

MM 



/Hv{t)]\ 



(8) 
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In the present case, the product rule for the functional derivative leads to 

ii^) c u\\ _ / ( 9i{x) \ 5x{t) \ 



>7»(a;) 

The response function Sx{t)/S£_i{t) is calculated from Eq. ([^ as \li 



(9) 



dx{t) ^ l g^{x) 

Sm 2 7.(:r)' ^ ^ 

where the + sign applies to i = 1 and the — sign to i = 2 since the relative coordinate has been defined as 
X = xi — X2 — xq. The factor 1/2 is due to the Stratonovich interpretation. The final result for the average over the 
multiplicative noise terms is 

1 //5^(^MM_M^)!2iM')\ (11) 



2\7i(2;) / 4\V 7^(2;)^ 
which can be substituted into Eq. The average center of mass velocity is then 



2\V7i(a;) l2{x) J I 4\ 7i(x)3 / 4^ ^^(a;) 



(12) 



Substitution of the g^, Eq. ([3]), yields our first main result, an exact expression for the average inchworm velocity in 
the overdamped regime: 



2\\7i(a;) l2{x) J I 2 \\^i\(x) 72(2;), 

A2 h2ix)'\ Ai /7i(a;)' 



4 \72(a;)V 4 \7i(a;) 



(13) 




The averages can be performed if the distribution of the relative coordinate x is known. Since the equation of 
motion for a;, Eq. ©J is decoupled from the center of mass coordinate, this distribution is determined by solving the 
Fokkcr-Planck equation corresponding to the Langevin Eq. ([6]) . In order to simplify notation let us rewrite Eq. ([6]) as 

x = a{x) + h{x)i{t), (14) 

where (^{t) is the superposition of ^i(t) and S,2{t) with the same Gaussian white statistics and, upon substitution of 
the gi, the auxiliary functions a{x) and b{x) are defined as 

(15) 
(16) 

The Fokker-Planck equation for the probabihty distribution p{x,t) associated with Eq. (fT4|) reads in Stratonovich 
interpretation [l9l | 

d d 1 d d 

-Q^P{x,t) = -—a{x)p{x,t) + 2 

If the averages in Eq. are calculated with p{x, t) as solution of Eq. p7|) . one obtains the full time-dependence of 
the average inchworm velocity. However, apart from the fact that the time-dependent solution is usually very difficult 
to find, we are also mainly interested in the steady state properties of the inchworm, where the center of mass is 
expected to move with a constant velocity. The solution in the stationary regime is determined in a straightforward 
way by setting dtp{x,t) — 0. Additionally, the probability flux can be set to zero since the relative coordinate is 
bounded in the spring potential. This simplifies the Fokker-Planck equation to 
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The solution is then found by straightforward integration and reads 

P(-)-A^^exp{/^^d,}, (19) 

where N is the normahzation constant. Eq. is our second main result, an exact expression for the stationary 
probability distribution of the relative coordinate. With this distribution we can explicitly calculate the average 
velocity of the inchworm in the steady state given by Eq. (fT3| . 



IV. COMPARISON WITH SIMULATION DATA 



In this section we compare the exact solution for the average inchworm velocity, Eq. (|13p . with the results from a 
direct simulation of the inchworm. In the simulation the coupled Langevin equations (|^|retaining inertia have been 
numerically time-stepped in the overdamped regime using an Euler-Maruyama scheme [20] . An asymmetric damping 
is assumed, where for simplicity the friction acting on bead 2 is set constant and the friction on bead 1 is chosen in 
the form of (see [lol |) 

7i(a;) = 1 + 7o w tanh(.T/w), (20) 

which prescribes a friction varying between two extremal values. The difference between these values is given by 2w"fo, 
where w parametrizes the width of the crossover regime and the sign of 70 specifies the orientation. For positive 70 
the damping increases for larger x and for negative 70 it decreases. 

For the parabolic spring potential U{x) — ^kx^ and friction term Eq. ([^0]) the distribution p{x), Eq. (|19p . has a 
lengthy but exact analytical form. With this distribution the averages in Eq. (fT3)l have been calculated by numerical 
integration. Fig. [2] shows a comparison of the exact solution with the simulation results of Fig. 1 in Theory and 
simulation show excellent agreement such that error bars have been omitted. For I70I > 0.22 a systematic deviation 
can be observed which is due to the fact that for the given form of the friction term, Eq. pp)) with w = 4.0, 1/71(2;) 
becomes singular if [70 1 approaches 0.25. This leads to an increased error in the numerical integration. 

As discussed in [ifll an interesting feature in Fig. [2] is the occurrence of current reversals for a variation of the noise 
strengths on each bead. In the case Ai = A2 the inchworm moves due to the asymmetric friction in the direction of 
the bead with constant damping, in agreement with the propulsion mechanism explained in Sec. [ill The direction 
of movement is reversed if the noise strength on bead 2 is increased. In that case the asymmetry of the noise input 
rather than the asymmetry of the friction dominates the movement. This becomes clearer below when the case of 
equal stretch dependent friction on both beads is considered. 

For another set of parameter values the distribution of the relative coordinate has been determined from the 
simulation and agrees well with the exact overdamped result Eq. (|19p (see Fig.[3|). 



V. SPECIAL CASES OF THE EXACT SOLUTION 



Having obtained an exact solution for the average inchworm velocity, it is possible to verify some of the perturbative 
results in [IQ] . We focus on relevant limit cases of the solution and highlight the importance of the nonequilibrium 
noise as well as the asymmetry and the stretch dependent friction for the movement of the dimer. 

1. On physical grounds it is obvious that for equal and opposite forces on the two beads no net movement should 
result. From the exact solution it is clear that for 71 (x) = 72(2:) and Ai = A2 all terms in Eq. (|13p vanish 
pairwise and therefore {xcm) = 0. 

2. In the absence of nonequilibrium driving the distribution is required to assume the equilibrium Boltzmann form 
and due to the second law of thermodynamics the average inchworm velocity should be zero. This can be seen 
as follows. For Ai = A2 = the distribution of x, Eq. p9)) . reads 



^ ^ V2fcBT(7i(2;)-i+72(x)-i) ' [J \ ksT 2 ^y)-' + j2{y)-' J "j ^ ^ 

Both integrals in the exponent are easily performed. Noting that the second integral yields a logarithmic term 
which subsequently cancels with the prefactor, the Boltzmann distribution p{x) readily recovered. 
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FIG. 2: The theoretical prediction Eq. (|13p compared with the simulation data of Fig. 1 in p^. We define /S.A = A2 — Ai, 
where the noise strength Ai is fixed at 1.0 and A2 is varied for the different curves. Note that in our notation 70 is 71 of 
Parameter values: k — 0.05, fcsT — 0.01, w = 4.0. 



Simulation 
-Tlieory 





FIG. 3: The exact expression for p{x) Eq. (|19|l compared with the distribution sampled from simulation data for small but 
finite inertia. Parameter values: k = 0.1, T = 0, ^2 = 0, 71 = 0.4, w = 2.0. 



Using the equilibrium distribution the first term in Eq. can be rewritten as 



7i(a;) 72(2;)/ / 2 7_oo V7i(a;) j2{x) J dx 

1 1 
71 (.t) 72(2;) 




(22) 
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where partial integration has been used in the last step. Therefore the remaining two terms in the expression 
for the average velocity cancel and we see that without nonequilibrium driving as expected (icm) = 0- 

3. The crucial ingredient in the inchworm model for the rectification of the diffusive motion is the stretch dependent 
damping. This is evident from the exact solution. For 71 and 72 both independent of x the distribution of x is 
p{x) oc exp{— CC/(x)/fcBT}, with a constant C. The average center of mass velocity is then given as 

(iem) OC {U'{x)) CX / —p{x)dx = 0. (23) 

This confirms that stretch dependent damping is required for directed motion of the inchworm. However, the 
asymmetry in the system does not necessarily have to originate from an asymmetric damping as discussed in 
the following. 



A. Asymmetric noise strengths 

In order to further investigate the influence of asymmetric nonequilibrium noise, we consider the special case of 
equal stretch-dependent friction on both beads ('symmetric friction'), 71 (x) = 72(2;), but asymmetric noise strengths 
■A-i 7^ ^2- In this case the distribution of the relative coordinate can be written as 

p(.)^^^I±^=exp(- r4L^;(^)7i(^) + 2W(y),l. 
v/4fcBr7i(x)+Ai+A2 \ J AksTjiiy) + + A2 " j 

The second integral yields a logarithmic term which can be absorbed into the prefactor. The result is 

^^^^°^7i(^) + (Ai+A2)/(4fcsr)'^"Pl J fcsT7i(y) + (Ai+A2)/4^^r ^^^^ 
In turn, the average inchworm velocity is obtained directly from Eq. (|13p as: 

{±c^)^-A^^){A^-A,l (26) 



4 \ ii{xf 

which is the exact analogue of the perturbative result Eq. (3) in . This becomes evident when we expand 71 (a;) in 
powers of x and truncate after the zeroth order. 

Assuming that the average in Eq. performed with the distribution Eq. (|25p is generally non-zero, we therefore 
establish that asymmetric driving in combination with stretch-dependent damping generates unidirectional center of 
mass movement. If we consider a fixed total noise strength Ai + A2 and vary the strength on each bead, we see 
that the distribution p{x) remains unchanged while the average inchworm velocity is exactly linearly proportional to 
the difference in noise strengths. It is then immediately evident that a sign change of A2 — Ai leads to a reversal of 
the direction of motion. This is a universal result for symmetric stretch dependent friction and independent of the 
particular form of the spring potential. 

A difference in the active noise on the two heads should be realisable experimentally in artifical systems. If we 
generalize the "Janus beads" of [21], [2^ by connecting two catalyst-coated beads of different size by a flexible 
polymer, we should obtain a self-propelling dimcr with asymmetric active noise. 



B. FENE spring potential 

In Eq. (pS]) it is indicated that with a continuous increase of one of the Ai the velocity of the inchworm is also 
increased indeflnitely. By contrast, a real molecular motor cannot hydrolyse an arbitrary amount of ATP and is 
expected to reach a saturation limit. In this shortcoming of the model has been explained by the unbounded 
spring potential, which is able to absorb an infinite amount of energy. Since our exact result for the velocity is valid 
for arbitrary potentials U{x), an investigation of more realistic spring potentials, as for example a finitely extensible 
nonlinear elastic spring (FENE), is possible and presented in the following. 

The FENE spring force is usually written as 
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where k denotes the spring constant and Xm is the maximal extension of the spring. For small extensions x, the 
FENE spring behaves like a linear Hookean spring. The restoring force is rapidly enhanced for increased extension 
and becomes infinite in the limit x Xm thus modelling the more realistic scenario of a finite extension. The 
corresponding potential is 

For the FENE potential the distribution of relative coordinate and the average inchworm velocity can be determined 
as before, via Eqs. and where simply U{x) has to be specified by Eq. (pS)) . In order to obtain explicit 

expressions it is necessary to perform the integrals numerically. We use the same form of the stretch dependent 
damping, Eq. (|20|) . and focus on two asymmetric scenarios, namely (i) asymmetric friction and equal noise strengths 
and (ii) equal friction on both beads and asymmetric noise. 

Fig.HI^a) shows results for the average inchworm velocity in case (i). We can see that up to Xm = 10 the behaviour 
of the quadratic spring potential is reproduced. Constraining the extension of the spring further, noticeably reduces 
the inchworm velocity. On the other hand, if the noise strength on each bead is equally increased (see Fig. [H^b)), 
even for small maximal extension the velocity is considerably boosted. In the symmetric friction case (ii) a similar 
observation is made in Fig.[n{a), where the velocity is plotted as a function of A2. Here the inchworm moves faster for 
smaller Xm- In all cases considered the qualitative features of the average inchworm velocity are very similar to the 
quadratic spring potential. A saturation of the velocity is not observed. Fig.l^Kb) shows a plot of the distribution p(a;) 
in the symmetric friction case. For Xm < 10 the FENE spring potential leads to an abrupt decay of the distribution 
in the vicinity of the maximal extension. 




FIG. 4: Average inchworm velocity for the FENE spring in the asymmetric friction case as a function of 70. (a) Noise strengths 
Ai = A2 = 1.0 and four different values of Xm- (b) Fixed maximal extension at Xm = 3.0 and increasing symmetric noise 
strengths Ai = Ai. Parameter values: k = 0.05, w = 4.0. 



VI. SUMMARY 



We have investigated a Brownian inchworm model of a self-propelled particle consisting of an elastic dimer driven 
by thermal and nonequilibrium noise. As our main results we derived exact expressions for the distribution of the 
relative coordinate and the average inchworm velocity for nonequilibrium Gaussian noise. The crucial property of the 
model allowing for an exact solution is the decoupling of the equations of motion of the center of mass and relative 
coordinate in the overdamped regime. The distribution of the relative coordinate is obtained in a straightforward 
way by solving the corresponding Fokker-Planck equation under the condition of stationarity and zero probability 
flux. For the average velocity it is necessary to determine noise averages over functionals of the noise, which can be 
calculated with the help of Novikov's theorem for Gaussian processes. 

The exact solution shows excellent agreement with results from a direct simulation of the equations of motion 
retaining inertia and provides the exact foundation for some of the perturbative results in [l^l- For a FENE spring 
potential the behavior of the inchworm does not change significantly compared with a harmonic potential and in 
particular the velocity does not saturate. An inchworm- type model on a discrete lattice has been discussed in [24| . in 
which the mechanisms of walking and inchworming are compared and criteria offered to distinguish which of these is 
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FIG. 5: The inchworm with FENE spring potential in the symmetric friction case, (a) Average inchworm velocity as a function 
of A2, where Ai — 1.0 and 70 = 0.1. (b) Distribution of the relative coordinate for Ai = 1.0 and A2 — 2.0. Parameter values: 
K = 0.05, w = 4.0, 70 = 0.1. 



operating in an experimental system. The internal states in this model are discrete and finite in number. It would be 
interesting to know whether the model, unlike ours, shows saturation in the velocity as a function of an appropriately 
defined parameter corresponding to ATP concentration. 

Our Brownian inchworm exhibits a variety of different movement mechanisms depending on the stretch dependent 
friction, the noise strengths on each bead and possibly the spring potential and statistics of the noise. So far only 
the case of Gaussian nonequilibrium noise has been thoroughly examined. Future work should include noise with 
non-Gaussian statistics, for example Poissonian shot noise, which might be a more appropriate way to model chemical 
energy input. In a regime where the frequency of stochastic kicks is low compared with the dimer relaxation rate, 
the inchworm can fully relax in between noise inputs and might reveal a qualitative different behavior. However, the 
treatment would then require diff'erent analytical methods and is left for future work. 
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APPENDIX A: ADIABATIC ELIMINATION 



In this appendix we perform the adiabatic elimination of the velocity for position-dependent friction. This derivation 
closely follows [T3]. 

Consider the following Langevin equation for a Brownian particle 



X = V, 

mv = —U'{x)--j{x)v + g{x)C{t), 



(Al) 
(A2) 



with — 0; {({0)({t)) = 6{t) and the prime denotes differentiation w.r.t. x. The Fokker-Planck equation 

associated with the above Langevin equation is 



dP{x,v,t) _ ^dP{x,v,t) ^ 1 d 
dt dx m dv 



U'{x) +^{x)v + 



d [g{x)f 



dv 2m 



P{x,v,t). 



(A3) 



Defining 



k(,x,t) = / dvvT{x,v,t) 



(A4) 
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we need to find the evolution equation for Qo{x,t). It is easily shown that 



dQo 


dQi 




dt 


dx 




dQi 
dt 


dQ2 
dx 


m 


dQ2 
dt 


dQ3 

dx 


27(2;)^ 

W2 

m 



U'ix) 



The above equations are exact. We now take the ovcrdampcd hmit (i.e., m/j{x) d/dt) 

m dQ2 U'{x) 



Qi 
Q2 



7(a;) dx 7(0;) 
m dQ3 U'{x 



2'y{x) dx 7 (a; 
Using the second equation above in the first, we get 



Qi 



2m'y{x) 



Qo- 



= ^ 9f m dQ.^U'jx)^^ [^^"^^]' q,^ ^'^""^Qo 
^{x) dx\2'y{x) dx 7(0;) 2m'y{x) °J 7(3;) ° 



1 d [gjx)?^^ U'jx)^^ 



2j{x) dx 7(0:) 



7(0:) 



to lowest order in 772/7(0;). Thus 



dQo{x,t) ^ d fu'jx) ^ 1 d l9i^)?]Qj^A 
dt dx \ ^{x) 2'y{x) dx ^{x) I " ' 



^ d U'jx) y(.rj[ff(.r;)]^ 1 d 
dx y -yix) 2[7(a;)]3 2dx 

This implies the following Langevin equations 

U'{x) i{xMx)r , 9{x) 



\Qo{x,t). 



X = ^ ■ ■ ; - - - |- 

7(a;) 2[7(a;)]3 7(3;) 

U'{x) g{x)g'{x) , g{x) 



X = —- 



j{x) 2[j{xW 
where 77 is a zero mean Gaussian white noise. 



ltd 

Stratonovich 



(A5) 
(A6) 
(A7) 

(A8) 
(A9) 



(AlO) 



(All) 

(A12) 
(A13) 
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